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TURBULENT AND LAMINAR HEAT -TRANSFER MEASUREMENTS ON A 

1/6- SCALE NACA RM-10 MISSILE IN FREE FLIGHT TO A 

MACH NUMBER OF k .2 AND TO A WALL 
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By Robert 0. Piland, Katherine A. Collie, 
and William E* Stoney 


SUMMARY 


Aerodynamic-heating data were obtained between Mach numbers of 1.2 
and k .2 from a free-flight test of a l/ 6 -scale finless NACA RM-10. 
Transient skin temperatures were measured at six stations along the 
body. The maximum skin temperature recorded during flight was l400° R. 
The test of Reynolds number covered a range from 4.2 x I 06 to 27.0 x 10 ^ 

Temperature recovery factors were obtained for each station at a 
single time diiring the flight. The values agreed with Prandtl niomber 
to the one-third power within ±3 percent. 

The experimental turbialent Stanton numbers were in fair agreement 
with Van Driest 's theory (assuming transition at the tip) during the 
time the skin was heating. Ikuring the skin-cooling period that followed 
the experimental values were consistently higher than the theoretical. 

During the skin-heating period as turbulent heat-transfer coeffi- 
cients were being measured, a sudden "burst" of laminar flow occurred 
for 0.5 second. The laminar or transitional flow extended to 38 percent 
of the body length or a Reynolds number of 13*3 X I 06 . For a short 
time prior to, during, and after the "bxirst," the Mach number and 
heating conditions were such that the two-dimensional stability theorem 
would have predicted the possibility of an infinite length of laminar 
flow. At one other very short period of time while the model was in 
the stability region, there was no evidence of laminar flow. 

Laminar heat-transfer coefficients were measured again during the 
latter part of the test. Transition from turbulent to laminar flow 
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began at local Reynolds numbers of 5 to 4 x 10^. Poor agreement was 
obtained with laminar theory at the foremost station. The agreement 
became progressively better moving rearward on the body . 


nWRODUCTION 


Reference 1 presented drag data from a free-flight test of a 
1/6-scale, sting-mounted NACA RM-10 model. In assessing these data, 
it appeared as if transition might have occurred during the flight. 
Full-scale RM-10 heat-transfer models which allow assessment of the 
condition of the boundary layer (refs. 2 and 5) had previously been 
tested over this Mach number range (l.O to iv.O) but at much higher 
Reynolds numbers and no evidence of transition or laminar flow was noted. 
Because of the need for information on transition and heating data at 
all Reynolds numbers, it was decided to test a similar l/6-scale RM-10 
instr\amented to measure skin temperature at stations along the model; 
these temperatvire data, when reduced to heat-transfer coefficients, 
would indicate by their magnitude the nature of the boundary layer. 

The heating data obtained from this flight are therefore reported 
herein for a Mach number range of 1.2 to k.2 and a Reynolds number 
range of k.2 x 10^ to 27.0 X 10^. The maximum skin temperature recorded 
d\iring the fli^t was l400° R. The test was conducted at the Langley 
Pilotless Aircraft Research Station at Wallops Island, Va., in 
January 1955* 


SYMBOLS 


Q quantity of heat, Btu 

A area, sq ft 

h local aerodynamic heat-transfer coefficient, Btu/sec-sq ft-°F 

T temperature 

e emlssivity 

a. Stefan-Boltzman constant, O.it-855 X 10"^, Btu/ft^-sec-°R 

b 

k thermal conductivity of air, Btu-ft/sec-°F-sq ft 
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Tw 

Pw 

Ch 
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A,B 

R.F. 

Pr 

J 

g 

M 

R 

P 

t 

X 
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thermal conductivity of Inconel, Btu-ft/sec-°F-sq. ft 
thickness, ft 

specific weight of wall, Ih/ft^ 
specific heat of wall material, Btu/lh-°F 
Stanton number 

specific heat of air at constant pressiare, Btu/slug-°F 

density of air, slugs /cu ft 
velocity, ft/sec 
ratio of specific heats 
pressiare, Ib/sq. ft 

gas constant for air, I716 ft^/sec^-*^ 

constants in approximate equation for Cy 

recovery factor, . Q-X . T 
■^so 

Prandtl number 

mechanical eq.uivalent of heat 
standard acceleration of gravity 
I>fech number 
Reynolds niomber, 

d 

viscosity of air, slugs/ft-sec 
time from start of flight, sec 
axial distance along body from the nose, ft 

standard deviation (with subscripts to indicate quantity involved) 
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Oqjj maximum probable error in Cjj 

Subscripts: 
aw adiabatic wall 

o undisturbed free stream ahead of model 

so stagnation 

V just outside bovmdary layer 

w wall ( skin) 

s sting 


MODEL, INSTRUMENTATION, AND TEST 


The test model was a I/ 6 - scale, finless NACA RM-10, sting mounted 
on a carrier body which housed the instrumentation and sustainer rocket 
motor. A photograph of the test model is shown in figure 1. A sketch 
of the test model and the carrier body is presented on figure 2. The 
model skin was formed by spinning a l/52-inch-thick sheet of Inconel, 
then highly polished. The skin- thickness measxjrements were made on the 
finished model by using micrometers. The surface roughness was of the 
order of 10 microinches. There was a break in the model at station 15 
for purpose of assembly. Shortly behind the break there was a circum- 
ferential row of flat-head screws. 

Skin-temperature measurements were made at six stations along the 
body as shown in figure 2(b) by use of thermocouples. The thermocouple 
wire was no. 50 chromel-aliamel . The jiuictlon between the wires consis- 
ted of a bead of about 0.01 inch in diameter, formed by fusing the wires 
together and using the mercury -bath technique . Care was taken that the 
wires were not in contact except within the bead. The beads were fitted 
into holes drilled through the Inconel skin at the proper stations, with 
the thermocouple leads inside the nose. The holes were then welded 
closed with Inconel welding rod and the exterior s\rrface was smoothed 
and polished. The cold junctions of the thermocouples were potted in 
paraplex Inside of a brass block of sufficient mass that no change in 
cold-junction temperature would occur during the relatively short time 
of the test. The cold- junction temperature was measured just prior to 
launching by a resistance-type teirqperatiore pickup also potted inside 
the brass block. A cross check of this temperature was obtained by 
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simultaneously measiiring the skin temperature with a thermocouple taped 
to the exterior surface of the nose . 

During flight, three standard voltages and the outputs of the six 
thermocouples were commutated and transmitted on a single telemeter 
channel. The commutation rate and the electronic system were such that 
each thermocouple voltage was transmitted l 4 times per second, and each 
standard voltage was transmitted 7 times per second. The three standard 
voltages, supplied by a mercury cell and a voltage -divider network, were 
chosen equivalent to the lowest temperature, the midrange teiirperatiire , 
and the highest ten 5 )eratirre that the skin thermocouples were expected 
to reach. Commutation and transmission of these known voltages along 
with the voltage readings of the skin thermocouples provided an in- 
flight check calibration of the thermocouple telemeter and recording 
system. 

The sustainer rocket motor was a 6.25-inch ABL Deacon rocket motor. 
The booster consisted of a cluster of three of these rockets. The model 
and booster on the launcher are shown in figure 3 • 

The model was launched at an angle of 70° to the horizontal. The 
booster propelled the model to a Mach number of 2.2. After a 1-second 
coasting period, the sustainer motor fired and accelerated the model 
to a peak Mach number of A. 2. 

During the time of rocket-motor firing and the coasting period 
that followed, an MCA telemeter located in the carrier body was relaying 
the temperatxire measurements to ground-receiving stations. The model 
velocity was measiared by CW Doppler radar, and SCR tracking radar 
measured the flight path, giving horizontal range, altitude, and flight- 
path angle. Atmospheric conditions were measured by radiosondes, one 
launched some time prior to flight to reach high altitude by the time 
the model was fired, and one at the time the model was launched to 
accurately measirre low-altitude data. Although no instrumentation was 
carried in the model to make a direct measurement of angle of attack, 
the size of the stabilizing fins was such as to make the carrier body 
a highly stable vehicle so that an angle of attack near zero (± 1 °) was 
maintained. Measurements such as flight path or static pressure (not 
presented herein) gave no Indication that the model followed other 
than a zero- lift trajectory. 


DATA REDUCTION 


During the flight test of the model, the following information was 
obtained as a function of flight time: 
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(1) Skin temperature measurements (fig. 4 ) 

(2) Model velocity, Mach nvunber, and altitude (fig. 5 ) 

(3) Air properties at any given time (density, static temperature, 
and speed of sovind) 

It is desirable to reduce this information into the form of Stanton 
number and recovery factors . The method of reduction used is described 
in the following sections; 


Stanton Number 

The basic heat-transfer eq.\iations as given in reference 4 are; 
For convection; 


^ 


For radiation outward: 


dt 


- A 


For radiation inward (model skin to internal sting) : 

(t^^ - Ts^)a 



1 

dt ^ 

+ i - 1 


^w 


For conduction longitudinally through skin; 


^ ^ A 

dx2 


( 1 ) 


( 2 ) 


( 5 ) 


( 4 ) 


Equations (5) and ( 4 ) are strictly applicable for the case of two coaxial 
cylinders and a cylinder, respectively. Since the model shape in the 
present case is parabolic, the equations are approximations. The error 
introduced in the heat-transfer coefficients by these approximations is 
discussed in the section "Accuracy of Stanton Numbers." 
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The time rate of change of heat contained in the skin is 


dt 




d^ 

dt 


A 


(5) 


The summation of the quantities of heat contributed by convection, 
radiation, and conduction (eqs. (l), ( 2 ), ( 5 ), ( 4 )) will equal the heat 
contained in the skin (eq. (5))* Neglect conduction, which is less than 
1 percent of the total heat transfer, and cancel A throughout; then 



From the local convective heat- transfer coefficient (eq. ( 6 )), the 
Stanton number can be determined: 


Ch = 


"=p,vPvVv 


(7) 


The properties of the test material (inconel) are known. The 
thickness Tv was measured with micrometers at the various stations 
and was between 0.051 and 0.055 inch. The density is constant. 

The specific heat varies with skin temperature, as shown on fig- 

ure 6 . These data were measured by the National Bureau of Standards 
and are presented in reference 5 . The emlsslvlty of the material is 
considered to be 0.5 throughout the flight. Measurements in reference 5 
show that for nonoxidized Inconel, the emissivity varies only slightly 
from the value of 0.5 for the temperature range of the test. The rate 
of change of temperature with time dT/dt was mechanically read, using 
the measxired temperature data. The sting temperature Tg used in 
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calculating the radiation term was assumed to he constant and equal 

to ground air temperature. The thermodynamic properties of air 

(cp, Pr, etc.) used in the computations were obtained from reference 6 

and are presented in figure 7. The adiabatic wall temperature Taw 
was obtained from the definition of recovery factor (ref. 5): 


R.F. 


'^aw ~ 

Tso " '^v 


This leads to 

Taw = Tv + Pr ^ ^Tgo - Ty^ 

The turbulent recovery factor is taken equal to Pr^/^ and the laminar 

recovery factor is taken equal to Pr^/^^ both based on wall temperatirre . 

The tiorbulent recovery factor was assimed equal to Pr^/^ because Pr^/5 
has been shown to approximate the meas\rred recovery factor within one 
percent for a variety of test conditions, (See ref. 7-) The measured 
recovery factors were not used since they were obtained at only one 
time during the flight and hence at only one temperature condition. 

The laminar recovery factor was assumed equal to Pr^/^, which is also 
a reasonable approximation as shown by reference 7 • The stagnation 
temperature was computed from the energy equation 



dT 


which takes into account the variation of the specific heat of air with 
temperature. Figure 8 gives a typical time history of adiabatic, stag- 
nation, and skin temperatures. The temperature Tq is known from model 

altitude and radiosonde measurements. The local temperature just out- 
side the boundary layer Ty is obtained by correcting the free-stream 
temperature Tq for the effect of body pressures at the various sta- 
tions, The slender -body theory of Jones and Margolis (ref. 8) was used 
for calculating the pressure distribution at supersonic speeds. Local 
values of velocity and density were also obtained for use in reducing 
the local convective heat-transfer coefficients h to Stanton num- 
bers (Ch = ^ ^ ' 

\ ^p,v^vV 
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Recovery Factor 


Recovery factors were determined at each tengjeratxrre measuring 
station at a single time during the flight by using the following method. 

If the convective heat transfer is assumed to be zero, equation (6) 
may be rewritten 




^w'^b^w^ + % 

1 

m ^ T ^ 

dTw _ 


i + 1 


dt 



Pw'’’w‘^w 



By a cut-and-try process, the slope dT/dt which satisfied the conditions 
on the right-hand side of the equation was found. The time corresponding 
to this slope will be the time of zero convective heat transfer. Now at 
the time of zero convective heat transfer, the temperature of the wall 
Tw will equal the adiabatic wall temperature Taw Therefore, at this 
given time we may determine a recovery factor, since 


R.F. 


^aw ~ ^v 
^so " ^v 


Accuracy of Stanton Numbers 

The accuracy of the Stanton numbers is dependent upon the possible 
errors in the various measirred and estimated quantities used in equa- 
tions (6) and (7). An appendix is devoted to a discussion of these 
errors and their probable resulting errors in Cg. Before discussing 
these more specific so\rrces of error, let us consider the applicability 
of the equations used in determining Cg. Equation ( 5 ) is strictly 
applicable for two coaxial cylinders whereas in the present case the 
model is parabolic in shape. In order to estimate roughly the order of 
magnitude of error that might be expected from this approximation, the 
relative contributions of convective and radiative heat transfer were 
determined. Figure 9 presents the variation of h with time for a 
typical station indicating the relative magnitude of the convective and 
radiative heat transfer. It is seen that the total radiative heat trans- 
fer is at most 20 percent of the total. Furthermore, the radiative 
heat transfer inward, not shown on the figure, was at the most only 
1.5 percent of the total heat transfer. Consequently, even if the use 
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of equation ( 5 ) results in a fairly large percentage error in the cal- 
culation of the inward radiative heat transfer, the error contributed 
to h would still be small. 


Equation (4) is strictly applicable for the case of conduction 
longitudinally through a cylinder. The heat transfer for the present 
parabolic case calculated from equation (4) was a negligible part of 
the total heat transfer. In addition, the longitudinal conduction 
through a conical skin was calculated and this contribution to the total 
was also negligible . 


Using the method presented in the appendix, the errors in Cjj 


attributable to errors in the various meas\ared and estimated quantities 
of equations (6) and ( 7 ) have been calculated for station 6. The meas- 
ured value of Cg, together with an accuracy band, are presented in fig- 
ure 10 as a function of time. The ratio of the probable error in Cg 
to the experimentally determined Cg, plotted as (JCg^Cg, is also pre- 


sented in this figure. It can be seen that the accuracy of Cg varies 
greatly during the test and this is primarily due to the variation of 
Taw - Tw (See fig. 10.) The greater the magnitude of this quantity, 
the greater is the accuracy of the values of Cg; and conversely, as 
Taw “ Tv approaches zero, the inaccuracy of Cg becomes so great that 
the data are of no value. The data presented in this report, therefore, 
include only those which are accurate within 20 percent. The accu- 
racy band and ratio presented in figure 10 are strictly appli- 

cable to the data at station 6. However, at the other temperature 
measuring stations the accuracy of Cg at any given time is essentially 
the same as the accuracy of the data at station 6. 


RESULTS AND DISCUSSION 
Experimental Recovery Factors 


The recovery factors determined as described previously are pre- 
sented in figure 11 as a function of model station. It should be noted 
that the recovery factors were evaluated at different times at different 

stations. The theoretical values of recovery factor Pr^/^ and Pr^/^, 
based on the wall temperat\ire at the times the recovery factors were 
determined, are also presented. The experimental values are- seen to 

agree with PrV3 within ±5 percent for all stations . Consequently, 
a turbulent boimdary layer is indicated at all stations at the time for 
which the recovery factors were computed. 
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Stauiton Number 

The experimental and theoretical (Van Driest, refs. 9 and lO) 

Stanton nxmibers are presented in figure 12 as a function of flight time 
for six temperature measuring stations. The parameters Ry, My, and 
Tw/Ty, which theoretically determine Cjj, are also presented on figure 12 
as a function of time. For all six stations the theoretical variation 
of Ch with time is presented. These calculations are obtained by 
using conical theory and assuming turbulent flow to exist from the nose 
tip. The theoretical flat-plate values have been modified in accordance 
with references 11 and 12 to give the proper value for a conical shape. 
The test model, although not a cone, is more nearly approximated by the 
cone than by the flat plate, at least at the forward stations . In addi- 
tion, at the most forward station a theoretical c\irve is presented where 
conical theory is again used and transition is assumed to occur at a 
Reynolds number of 2 x 10^. At the most rearward station a theoretical 
curve is also presented where calculations are based on flat-plate 
theory . 

It should be kept in mind that the theories of Van Driest with 
which the experimental data are compared were developed on the assump- 
tion of a constant wall temperature and constant pressure distribution. 

As can be seen in figure 4 there is considerable variation of tempera- 
ture along the body at any given time. The pressures also along the 
body are not constant, the body being neither a flat plate nor a cone. 

The effect of the pressure gradient is taken into account in obtaining 
the Stanton numbers by using values of density and velocity just out- 
side the boundary layer. 

Consider first the data obtained while the skin is being heated. 
These data cover the period of accelerating flight and a small portion 
of decelerating flight and extend to a time about 10 seconds after take- 
off. Except for a short period of time between 5.5 and 6.0 seconds, the 
experimental Stanton numbers (fig. 12) are in fair agreement with the 
turbulent theory res\ilts obtained by assuming transition at the nose tip. 
From station to station the agreement varies somewhat, but no signifi- 
cant consistent trends away from the theory are noted. At station 6 
the theoretical values of Cjj obtained on the assvunptlon of transition 
at a Reynolds niimber of 2 x 10^ do not differ greatly from those obtained 
on the assumption of transition at the nose tip. Up to about 8 sec- 
onds after take-off the experimental data are in better agreement with 
the theory based on the assumption of transition at the nose tip . 

Between 8 and 10 seconds after take-off the agreement is better if a 
transition Reynolds number of 2 x 10^ is assumed. 

At 5*5 seconds after take-off at stations 6, 7*7, and. S-k, with 
the local Reynolds number varying from x 10^ to 6.5 X 10^, the Stanton 
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nuiribers dropped, "to the theoretical laminar level. This result must 
have been due to a sudden "burst" of laminar flow. The effect was most 
pronounced at station 6 and decreased at the more rearward stations . 

No effect can be seen at station 11.1. In an effort to determine the 
cause of the laminar burst, recourse was made to the two-dimensional 
laminar boundary-layer stability theory of reference 15 . Figure 15 
presents a region defined by values of temperature ratio Tw/Ty and 
Mach nvimber in which, theoretically, infinite stability of the laminar 
boundary layer is attainable. The flight data show the model to have 
been in this region twice, between Mach nimibers I.5 and 2.0 and between 
Mach numbers 2.5 and 5.8- It was during the second period in the sta- 
bility region that the burst occurred, as shown in this figure. 

Between about 10 and 12 seconds after take-off, the skin temperatures 
reach maximums and the forcing function T^^ " Passes throu^ zero. 

As discussed in the section entitled "Acc\aracy of Stanton Numbers,' this 
phenomenon resilLts in inaccuracies in Cg large enough to make data in 

this interval meaningless and they therefore are not presented. 

After about 12 seconds after take-off the model is decelerating, 
and Mach number, Reynolds number, and wall temperatures are decreasing. 
The measured Stanton numbers are seen to be consistently and considerably 
higher than values predicted by turbulent theory (transition at nose 
tip) \mtil the time when transition to laminar flow begins (approx. 

15 sec) . For a possible explanation refer to figure 12(a) which pre- 
sents the data obtained at station 6. It can be seen that in this time 
interval the experimental data are in fair agreement with the theoret- 
ical values when transition is assimed at a Reynolds nvimber of 2 x 10 . 
Since transition occurs shortly afterward at Reynolds numbers of about 
5 X 108, this explanation of the seeming disagreement is likely. It 
should also be kept in mind that the accuracy of the data in this time 
interval is about ±12 percent which co\ild account for a part of the 
disagreement . 

6 

Transition is seen to begin at a Reynolds number of 5 to 4 x 10 
at about 15 seconds. Laminar heat-transfer coefficients are obtained 
as far rearward as station 11.1 (figs. 12(a) to 12(d)). Station 12. o 
(fig. 12(e)) shows the beginning of transition, but a laminar level 
had not been reached at the latest time (20 sec) for which data are 
available The laminar coefficients at station 6 are not in good 
agreement with the laminar theory. At the more rearward stations fair 
agreement is obtained. The possibility exists that the distribution 
of temperature along the nose, forward of station 6, is of such a non- 
uniform nature due to the particular construction of the model as to 
affect the measured laminar Stanton numbers . At more rearward stations 
on the body the nonuniform temperature near the nose tip would probably 
have less effect. 
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Figure 14 presents the variation of Stanton number with station in 
inches from the nose tip at various times dinring the flight. At the 
time of 1.8 seconds, the measured Stanton numbers are lower than the 
turbulent theory. Possible inaccuracies in the data at this time could 
account for most of the disagreement. At the time of 4.2 seconds the 
agreement with theory is good. At the time of 5*8 seconds the typical 
laminar, transition, turbulent curve is obtained. This time corresponds 
to the time of measuring the lowest heat-transfer coefficients during 
the laminar burst. At 7 .O and 9*0 seconds, the data are again in good 
agreement with turbulent theory. At l4.0 seconds, during skin cooling, 
the poor agreement discussed previously is obtained. 

Assuming a Reynolds number of transition of 2 x 10^ improves the 
agreement, especially at station 6 , as has been previously discussed. 

At the more rearward stations, however, the agreement is not iniproved 
significantly. At 18 seconds the variation of Cji with station indi- 
cates the laminar, treinsition, and turbulent regions . The laminar coeffi- 
cients, however, are considerably higher than the theoretical laminar 
values and the turbulent value at station I 8 is somewhat higher than 
turbulent theory. At 20 seconds the laminar coefficients, with the 
exception of station 6 , are in fair agreement with the theory. As stated 
earlier, in considering the laminar coefficients it should be remembered 
that the theory used here was developed on the assumption of a constant 
temperature ahead of the measuring station. 


CONCLUDING REMARKS 


An NACA RM-10 l/6-scale model has been flight-tested and convective 
heat-transfer coefficients (Stanton numbers) have been obtained at six 
stations along the body. The Mach number range was from 1.2 to 4.2 and 
the Reynolds number range was from 4.2 x 10^ to 27.5 x 10^* The maximum 
recorded temperature was l400° R. The following observations have been 
made from a congiarison of experiment and theory: 

1. The Stanton numbers predicted by the turbulent theory of Van 
Driest (on the assumption of transition at the nose tip) were in good 
agreement with experiment during the period when the skin was being 
heated. During the cooling period that followed, the measured Stanton 
numbers were consistently higher than theoiy. 

2 . The measured recovery factors obtained at each station at a 
single time during the flight agreed with Prandtl number to the l/3 
power within ±5 percent, indicating that turbulent flow existed at the 
time at which recovery factors were determined. 
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3 . A sudden "burst" of laminar flow occurred at a Mach number of 
2.9 and extended to 38 percent of the body length corresponding to a 
Reynolds number of 13-3 X lO^. The model at this time was in the theo- 
retical region of infinite stability of the laminar boundary layer. 

4 . Transition from turbulent flow begins toward the end of the 
test when the Reynolds number has dropped to 3 to 4 x I 06 . At this time 
the model was well out of the region of laminar boundary-layer stability. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February I 6 , 195^. 
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APPENDIX A 

ESTIMATED EE®ORS 
By William E. Stoney 


Accxaracy of Cg 

Errors may be grouped under two headings, systematic and random. 
(This discussion follows the treatment in ref. 14 .) Random errors are 
those which in a large number of measurements are as often negative as 
positive, and thus they affect the arithmetic mean but little. All 
other errors are classed as systematic. 

Systematic errors can be evaluated only by comparison of the final 
data with previously substantiated results, either theoretical or empir- 
ical. Although the number of previous tests with which such congjarlsons 
can be made is not large, such comparisons have indicated that there has 
been little or no systematic error present and such will be assiamed in 
the present case. 

Random errors appear in all the measured quantities. These indi- 
vidual errors will be quoted as standard deviations (cr), that is, the 
root mean square of the deviations of a set of observations from its 
mean value. In general, these values will be estimated from what are, 
statistically speaking, small samples and are thus themselves approxi- 
mate in nature. The effect of the individual values of a on Cg will 
be presented in the form 


Cg dx Cg 


As the Individual errors act independently of one another, their combined 
effect on Cg will be given by the following equation: 



(Al) 
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The equation for Cg reduced to the quantities which are measured, 
calculated, or assmed is presented below: 



where the approximate expression (A + BTv) has been substituted for 
c^. In the actual data reduction, experimental values of c^ were 


used. 


Also, the approximation 


1 + 




has been substituted for 


the equation 


2Jg 



dT 


which was used to get Tgg in the text. Equation (A2) will be referred 
to ocasionally in its more compact form: 


C = PwV^w ^ 

Pv^v*^p, v(^aw “ 

where r refers to the radiation terms. 

The equations for due to the errors in the individml quan- 

Ch 

tities are presented below. The e3<q)ressions for those quantities which 
have been assumed or which have negligible effect on the acc\aracy of 
Cg have been omitted. 
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Queintlty 


due to a 
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/ 50 \PTv/To 

(Taw “ 

■‘■0 


V^aw - V^v/To 


“W 



cr c T ^ 

T ^ - T \ 




Pb^w^w 

")/ 

dev 



r 




Pew 

n *T r* - 

dT^ 



ew 


dt 

+ 1 



r 





When values of the Individual errors (see table at end of the 
appendix) are substituted in the approximate forms of the error terms 
and the results s\jmmed as in equation (Al), the following equation is 
obtained ( for Mq = ^ ) : 


*^^Htotal 


^03 

‘ (Taw - T.w )2 



0.0022 


1/2 


As shown by this equation, the data become extremely inaccurate when 
dT 

Taw - Tv and ^0; this effect is compounded since both these terms 

dt 
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go to zero at about the same time. A comparison of this equation with 
equation (Al) computed using the complete error terms is shown in fig- 
ure 10. 

As can be seen from the upper part of figure 10, the percentage 
error where the difference is over 100^ R is of the order of 

20 percent or less. The percentage errors in this upper plot have been 

used to calculate values of Octj which have been added and sub- 

^Htotal 

tracted from the data and the results presented as a shaded area in the 
lower plot. It seems apparent from the comparison of the shaded band 
and the theoretical curve that nowhere are the data in significant dis- 
agreement with the theoretical calc\Hations . This is more apparent 
when it is realized that the values from theoretical calculations should 
be even higher than shown here due to the presence of laminar regions 
on the body. (See fig. 12(a).) This conclusion applies to the data 
of the remaining stations as well. 


Estimation of Errors of Individual Measurements 

Wall and shield temperatures .- Assuming the absence of systematic 
errors in the temperature measurements permits the random errors to be 
divided into two classes. In one group, the errors are such that they 
do not vary rapidly with time - for example, calibration errors and 
slow telemeter drifts. The other type exhibits itself in the data as 
random scatter with time and is probably due in great part to the errors 
involved in the recording and reading of the data. 

This latter type of error can be evaluated from the data when 
dT/dt is less than about 25^ per second. Such evaluation gives aj = 5^ 
This error is important only in flight conditions which give high rates 
of temperature increase but even there it is overshadowed by the first 
type of error, as will be shown. 

The error of the first kind can be evaluated only by comparisons 
of tests made imder identical conditions with different telemeters . 

Since this comparison has not been made (in fact it would prove very 
difficult to make in the temperature ranges of interest because of the 
difficulty of knowing or reproducing the conditions), it is necessary 
for a value to be assumed on the basis of experience with past measure- 
ments. In general, the accuracy of telemeter data of all kinds has 
been roughly estimated as 1 to 2 percent of the full-scale range of the 
instrument used. Previous tests using temperature equipment similar to 
that used on the present tests have indicated that teir 5 )erature data are 
at least as good and are probably better than the larger of these figures 
Thus for the purposes of this evaluation a figure of 1 percent of the 
full-scale value will be used and this estimate gives a value of aip 
of 2 Q 0 . 
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Error in Pq and Tq.- Errors involved in the determination of 
both Pq and Tq are a combination of three independent factors: 

(1) Error in the measiared altitude of the model 

(2) Error in the calculated altitude of the radiosonde 

(3) Error in the measurement of the quantities Pq and Tq 

In the altitude range of the present tests the estimated accuracy 
of the two altitude measurements together with the rates of change of Pq 
and Tq with altitude give errors in these quantities of dp^ = O.OO6 

and dip =0.2°. Unpublished NACA tests give values of dp^ = O.OO83 

and =0.9° for the radiosonde instrument error in the altitude 

range below 50,000 feet. Since these errors are independent, they can 

be combined by the equation o^otal ~\f^a °b^' gives the total 

d as follows: '' 

dr, = 0.015pr, dm = 1° F 

Mo O 


Error in Vq.- In the altitude range considered the major error in 

velocity measurement is due to the inaccuracies in the measurement of 

wind velocity. Since the wind velocity is not measured at the time and 

place of the model, any accxaracy which may be quoted is open to some 

doubt. However, consideration of accviracy of the actual measurements 

used to obtain the wind velocity as it affects the radiosonde balloon 

leads to a figure of Oy = ^ feet per second. 

*0 

dT 

Error in — If the meas\ired tenqoeratures are assvuned correct, 
dt 

dT 

errors in — — are caused by inaccuracies in the fairing of the tempera- 
dt 

ture data and in the reading of the slopes of the faired cvirves . An 
inspection of the scatter of such readings leads to a figure of 

= 1° per second for the faired results. This value seems to 

dT 

include "both types of error since an average error of per 

dt 

second was obtained by integrating the slope curves and comparing the 
resulting temperatures with the original temperature curve for the data 
of a typical station. 


^dTy j dt 
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p V T 

Errors in and Errors In the quantities pip, 

^o ^o _o ° 

and T^/Tq may all be referred to error in the calculation of 

the pressures over the body. While this error will be a function of 
the computation methods used, it can be seen from the comparison of the 
various methods with the results of the method of characteristics in 
reference 15 that a value of cTcp = 0.05Cp is not unreasonably low. 

IThis comparison of course includes the assumption that the method of 
characteristics gives correct values, and this assvnnption has been shown 
to be a good one except for fli^t conditions where the boundary layer 
mi^t be expected to appreciably affect the flow contour about the body. 
This latter condition should not exist in the present tests. Using the 

TT /tt fTi cant be obtauxed as 

o^''v/''o^ v/'^o 

functions of P^^'Pq* This value is, for any particular model, a function 

of Mach number and body station. However, for the present body and Mach 
numbers, Py/Po is small and leads to values of 


value Cq = 0.05cp values of 


rl^ 



^ 0.01 

Po 



T 

±SL 


Error in t^.- The accuracy of the measurement of the wall thickness 

is estimated to be about 0.001 inch. For the wall thickness of the 

present model this estimate leads to a value of a_ « 0.05t„. 

I w ” 

Error in €y. - The error in this quantity was estimated from the 

apparent scatter in the data of its measured values and is ct_ =0.02. 

^w 

Error in Cp^, p^, ct^, c^., and R.F..- The values of Cp^, 

'^b^ °w> R.F. have been assumed to be constants in which no enror 

appears. Although it is assumed that c^ is known absolutely, it is a 
quantity which varies with T^ and the error attributable to this fact 
appears in the error for T^. While Cp^ is a function of T^, T^ 
varies so little that any error due to an error in Ty is quite 
negligible . 

In the presentation of the data the recovery factor R.F. is 

assxmied to be equal to and Erl/2 for turbulent and laminar 

boundary layers, respectively, and thi:is the errors are assumed 0. It 
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shoiild be noted, however, that the equation for the error in Cg due 
to R.F. error, 


^ ^av ~ °R.F. 

\CH /R.F. Taw - Tw R-F- 

shows that small errors in this quantity can lead to quite large errors 
in Ch when 0. 

The values of the errors used in obtaining the data of figure 10 
are presented in tabiilar form as follows: 


Quantity 

a 


20^ R 

Ts 

20° R 

Po 

0.013Po 

To 

H 

0 


A fps 

dTw/dt 

1 ° per sec 

■^w 

0.05t^ 

p„ T„ V„ 

P,r T V„ 

V V 

^ 0.01 — . — , 

V ^ 

Po 0 ''0 


0.02 
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Figure 1.- Photograph of l/6 scale, sting-mo\mted NACA RM-10 test model. 


ro 

VJl 


NACA RM L 56 CO 5 



2J0 


rv) 

o^ 


(a) Model mounted on carrier body. 



(b) Model alone. 


Figure 2 . - Sketches of 
showing teraperatiire 


test model and carrier vehicle and of model alone 
pick-up stations. (All dimensions are in inches. ) 


NACA RM L56CO5 


NACA RM L56CO5 


27 



L-88158 

Figure J*- Model with booster on launcher ready for firing. 
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Figure 4 ,- Skin temperatures meas-ured for six .stations 
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Figure 5.- Measured flight conditions. 
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Figure 6 .- Instantaneous specific heat of Inconel. 
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Figure 7.- Thermal properties of air, reference 6. 
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Figure 8 .- typical time history of stagnation, adiabatic, and skin temperatures. 
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Figure 9 .- Heat- transfer coefficients for station 7.7 inches showing the 
effects of radiation and convection. 
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Figure 10.- The ratio of Stanton nxmber probable error to experimental 
Stanton number, heating potential and accuracy band for experimental 
Stanton ntunbers varying with time at station 6. 




Figure 11.- Variation of measiired and theoretical recovery factors with 

station. Prl/3 and Prl/2 are plotted for the times used in 
obtaining the measiared recovery factors. 
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Figure 12 .- Stanton nmber Cg and the quantities affecting Stanton number 

for the six temperature stations. 
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Figure 12 . - Continued 
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Figure 12.- Continued 
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(e) Station 12.8 inches. 


Figure 12.- Continued 
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(f) Station I8.O inches. 


Figure 12 .- Concluded 


■r 

ro 



Figure I5.- Con5)arison of flight conditions with conditions for infinite stability of the laminar 

boundary layer. 
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(a) 1.8 seconds. 


(b) 4.2 seconds. 




(c) 5-8 seconds. 


(d) 7-0 seconds. 


Figure l4.- Stanton number varying with stations along the body for a given time. 
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(e) 9«0 seconds. 


(f) l 4.0 seconds. 




(g) 18.0 seconds. 


(h) 20.0 seconds. 


Figiare l 4 .- Concluded. 
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